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RESEARCH MEMORANDUM 


A THECREIICAL INVESTIGATION OF THE DYNAMIC LATERAL 
OSCILLATORY STABILITY OF AN AIRPLANE 
HAVING A 60° TRIANGULAR WING 
By Joseph. L. Johnson 

SUMMARY 


A theoretical study has heen made of the dynamic lateral stability 
characteristics of an airplane having a 60° triangular wing. The calcu- 
lations included the determination of the neutral— lateral— oscillatory— 
stability boundary (R = 0), the period and time to damp to one— half 
amplitude of the lateral oscillation, and the time to damp to one— half 
amplitude for the spiral mode. Factors varied in the investigation 
were lift coefficient and altitude. Since small changes in some of the 
mass and aerodynamic characteristics may cause considerable differences 
in the lateral stability of an airplane, it is well to keep in mind that 
the results presented herein apply only to an airplane having the 
characteristics for which the calculations were made. 

The results of the investigation showed that the lateral oscillation 
of the airplane was stable over the lift— coefficient range at sea level 
and at an altitude of 30 , 000 feet. The damping of the lateral oscil- 
lation met the U. S. Air Force requirements for all conditions except 
In the low— lift— coefficient range at 30,000 feet. The airplane was 
spirally stable for all conditions. 


INTRODUCTION 


A theoretical study has been made of the dynamic lateral stability 
characteristics of an airplane having a 60° triangular wing. A three- 
view sketch of the airplane is shown in figure 1. The analysis was 
made by the Langley free— flight— tunnel staff and the calculations were 
made at the Langley Laboratory on a relay computer. 
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Calculations were made to determine the neutral— lateral— oscillatory— 
stability boundary (E = 0), the period and. time to damp to one-half 
amplitude of the lateral oscillation, and the t ime to damp to one— half 
amplitude for the spiral mode for the airplane over the lift— coefficient 
range with a gross weight of lh,llh pounds at sea level and at an altitude 
of 30,000 feet. The results of the investigation are presented in the 
form of stability charts where the boundaries E = 0 are plotted as 
functions of the effective dihedral parameter — ( and the directional- 

stability parameter C n g % The period and time to damp to one-half 

amplitude are presented as functions of lift coefficient. 


Recent theoretical investigations (reference l) have shown that small 
changes in some of the mass and aerodynamic characteristics may cause 
considerable differences in the lateral stability of an airplane and it is 
well to keep in mind, that the results presented herein apply only to an 
airplane having the characteristics for which the calculations were made. 


SYMBOIS AND COEFFICIENTS 


S 

■S' 

V 

b 

q. 

p 

w 

g 

m 

^b 

a 


wing area, square feet 
mean aerodynamic chord, feet 
airspeed, feet per second 
wing span, feet 

dynamic pressure, pounds per square foot 
air density, slugs per cubic foot" 
weight, pounds 

acceleration of gravity, feet per second per second 
mass, slugs (W/g) 

relative density factor based on wing span (m/pSb) 

angle of attack of reference axis (fig. 2), degrees 

angle of attack of-principal longitudinal axis of airplane , 
positive when principal axis is above flight path at the 
nose (fig. 2), degrees 


n 


NACA EM L9A014- 


3 


t 

P 

0 

E 


angle between reference axle and principal axis, positive 
when reference axis is above principal axis at tbe nose 
(fig. 2), degrees 

angle between reference axis and horizontal axis, positive 
when reference axis is above horizontal axis at the nose 
(fig. 2 ), degrees 

angle of flight to horizontal, axis, positive in a climb 
(fig. 2 ), degrees 

angle of yaw, degrees or radians 

angle of sideslip, degrees or radians 

angle of bank, radians 

Routhis discriminant (E* = BCD — AD 2 — B^E where A, B, C, 
D, and E are constants representing coefficients of the 
lateral— stability equation) 

radius of gyration about principal longitudinal axis, feet 




radius of gyration about principal vertical axis, feet 


% 




nondimens ional radius of gyration about principal longitudinal 

axls (*XoA ) 

nondimens ional radius of gyration about principal vertical 

aJClS (%fj 

nondimens ional radius of gyration about longitudinal stability 


axis ^\ykx o 2 cos 2 Tj + Kz o 2 sin 2 Ti^ 




nondimens ional radius of gyration about vertical stability 


axis 


(^f~ 2coe2, n + % 0 2sin2 ri.) 


tKZ 


nondimens ional product— of— inertia parameter 
I cos q sin 


((^o 2 -\^) COB 71 Binr] ) 
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lift coefficient (Lift/qS) 

yawing-moment coefficient" (Yawing moment /qSb ) 
rolling-moment coefficient (Rolling moment / q£Tb ) 
lateral— force coefficient (Lateral force/qS) 


rate of change of lateral— force coefficient with angle of 

sideslip, per degree or per radian, as specified (dCy/dp) 


rate of change of yawing-moment coefficient” with angle of 
sideslip, per degree or per radian, as specified (SC^Sp) 

rate of change of rolling-moment coefficient with angle of 
sideslip, per degree^ or per radian, as specified (dC^/dp) 


rate of change of lateral— force coefficient with rolling- 
angular^-velocity factor, per radian 


angular-velocity factor, per radian 



rate of change of rolling-moment coefficient with rolling- 

2vy 


rate-of change of yawing-moment coefficient with rolling- 
angular-velocity factor, per radian 



rate of change of rolling-moment coefficient with yawing— 

27/ 


angular-velocity factor, per radian 


rate of change of yawing-moment coefficient with yawing— 


angular-velocity factor, per. radian 
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Cy 

Z 
z 

P 
r 
3>b 
sb 
X 

t time, seconds 

P period of oscillation, seconds 

T-, /p time for amplitude of oscillation to change by factor of 2 

' (positive value indicates a decrease to half amplitude, 

negative value indicates an increase to double amplitude) 

C -^/2 cycles for amplitude of oscillation to change by factor of 2 

' (positive value indicates a decrease to half amplitude, 

negative value indicates an increase to double amplitude) 


rate of change of lateral— force coefficient with yawing- 
angular-velocity factor, per radian 

tail length (distance from center of gravity to rudder hinge 
line), feet 

height of center of pressure of vertical tail above fuselage 
axis, feet 

rolling angular velocity, radians per second 
yawing angular velocity, radians per second 
differential operator (d/ds-b) 
distance along flight path, spans (Yt/b) 
complex root of stability equation (c ± id) 



EQUATIONS OF MOTION 


The nondimsnsional lateral equations of motion (reference l), referred 
to a stability— exes system (fig. 3)j are: 

In roll 


2db^&x 2 I>b 2 0 + K XZPb 2 ^ = 


<V 


2 r 


+ 7pl r Vb* 


+ 
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In yaw 

- %P + *Pnjt0 * 


In sideslip 


2n b (DbP + T^t) = Cy p P + IpYjPbP + + iCy^t + (C L tan y* 


When is substituted for 0, ^ Q e^ a t for ♦ , and p 0 e*- s b 

for 0 in the equations written in determinant form, A must he a root 
of the stability equation 

AxA + BA 3 + CA 2 + DA.+ E = 0 (l) 


where . : ' 

A = % b 3 (Ky 2 K Z 2 - Kxz 2 ) 

B = -2M.-b 2 (W^E^C + % 2 C nr + - SKjz^Y^ ~ X JZ C l T ~ %Z C n p ) 

C = Hh^Ky^r^Cyp + ^u^Ky 2 ^ + Kz 2 C Zp Cy p + - ^TZfll^pY^ - W^XZCZp 

- CjipKjzCYp - ^CnpQXj, + KxzC np C Yp - Kz 2 Cy p C Zp - Kx 2 Cy r C np + KyzCy^j^) 

D ■“^ C n ar C Ip C Yp ” fA B C 2p C np + J C n p C l r C Yp + ^b c n p C Zp + 2 ^b c L%Z C ng 

- 2n- b 0 r ft 2 C Zp - ^^ Z 2 C np C L tan y + 2^^^^ tan 7 + jj=C z p C np C ^ 

“ iS 0 *p°rp - ^ r % C Y p + 

E = 2 ° L ( Cn T% ~ °lr Cn fi) + tan r^ZpCnp - 
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The damping and period of the lateral oscillation in seconds are given, 

^ p_2*£ 

c y a v 

where c and d are the real and imaginary parts of the complex root of the 
stability equation. The damping of the spiral mode is determined similarly 
from one of the two real roots (usually the less stable one) of the stability 
equation. 

The conditions for neutral— oscillatory stability as shown in refer- 
ence 2 are that the coefficients of the stability equation satisfy South's 
discriminant set equal to zero 


respectively, by the equations T, 


1/2 


E = BCD - AD 2 - B 2 E = 0 


and that the coefficients B and D have the same sign, 
sign of the coefficient B is determined by the factors 


In general, 
“ C Yp 3 ~ c n r > 


and 



which appear in the predominant terms of 


B. 


Thus, B is 


the 


positive in the usual case in which there is positive weathercock stability 
and positive damping in roll. Hence, the coefficient D must be positive 
if E = 0 isa neutral— oscillatory— stability boundary. 


Calculations 


Calculations were made to determine the neutral— lateral-oscillatory— 
stability boundary (E = 0), the period and time to damp to one-half 
amplitude for the lateral oscillation, and the time to damp to one-half 
amplitude for the spiral, mode over a lift-coefficient range (0.1, 0.5* 
and 1.00) for the airplane with a gross weight of 14,114 pounds at sea 
level and at an aiitude of 30*000 feet. 


The aerodynamic and mass characteristics used in the calculations 
are presented in table I. Yalues of C n and Cy R 

P (tail off) P (tail off) 

were taken from force tests made on a s i milar configuration in the 
Ames 40-hy 80— foot tunnel. The tail-off values of C^, C^ , C^, 

and Cj were estimated from reference 3 and from unpublished experi- 
mental data obtained in the Langley stability tunnel. The contributions 
of the tail to the stability derivatives were estimated from the equa- 
tions given in the footnote of table I and are s im ilar to those given in 
reference 4. The estimated values of C-n for the complete model for a 

lift coefficient of 0.1 and 0.5 are in reasonable agreement with 
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experimental data from the Langley st ability t unne l , hut no reliable 
experimental data are available for comparison with the calculated value 
for a lift., coefficient of 1.0. 

The roots of the stability equation were computed to. determine the 
period and time, to damp to ‘one— half amplitude for each condition for which 
a boundary was calculated. Presented in table II are the aerodynamic and 
mass characteristics used in computing the roots of the stability equa- 
tions. The values are essentially the same as those given in table I 
except that the mass characteristics are presented in a different form 
and the total derivatives are given for each condition. 


RESULTS AND DISCUSSION 


The results of-the calculations are presented in figures 4-to 10. 
All results are presented- in terms of lift - coefficient and figure . 9 , is 
presented for convenience in interpreting the results in terms of air- 
speed. Figure 10 is presented to show how the damping of the airplane 
compares with the damping requirements of the TJ. S. Air Force. 


Neut ral— Lateral— Osc illatory— Stab ility Boundaries 


The effect of lift coefficient on the neutral— lateral— oscillatory- 
stability boundaries (R = 0) Is presented in figures 4 and 5 for the 
airplane at _ sea level and 30., 000 feet, respectively. The symbols on 
these figures represent the C n and —C values for the airplane at 


P 


each lift coefficient. Figure 4 shows that _ the airplane was oscillatorlly 
stable for each lift coefficient and that the boundary shifted downward 
with an increase in lift coefficient - ; 


Increasing the altitude to 30,000 feet produced very little change 
In the position of .the .JR = 0 boundary for a given lift coefficient as 
shown by a comparison of figures 4 and 5» A small destabilizing shift 
of the boundary "occurred atr - = "0.10 and was the largest change in 
the boundary over the lift— coefficient range'. 


Period and Damping of the Lateral Oscillation 

The effect of lift coefficient on the reciprocal of the time to 
damp to oner-half amplitude and the period of the lateral oscillation is 
presented in figures 6 and J, respectively, and the results are tabulated 
In table HI. The airplane Is stable for all conditions of altitude and 
lift coefficient with the greatest damping occurring at sea level. The 
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results of figure 7 show very little change in the period of the 
oscillations with change in altitude. 

Presented in figure 8 are the damping characteristics of the air- 
plane in terms of cycles. The cycles to damp are determined by dividing 

the time to damp by the period. If the value of — - — is greater 

c l/2 

than 0.50 it means that the lateral oscillation damps in 2 cycles or 
less. It is seem from this figure that the lateral oscillation damps 
in 2 cycles or less for all conditions except at a lift coefficient of 
0.10 at an altitude of 30,000 feet. 

Presented in figure 10 are the calculated damping characteristics 
for the airplane in comparison with the required damping specifications 
of the TJ. S. Air force (reference 5). The symbols in this figure 
represent the points under consideration for the airplane and show that 
the conditions investigated meet the Air force requirements for satis- 
factory damping except at the low lift coefficients at an altitude of 
30,000 feet. 


Spiral Stability 

The reciprocal of the time to damp to one-half amplitude for the 
.spiral (aperiodic) mode is presented in table TTT , It is seen that the 
spiral mode is slightly stable over the low and medium lift— coefficient 
range. At a lift coefficient of 1.00 the aperiodic mode disappears and 
a stable long-period oscillation occurs. The sea— level condition showed 
the greatest amo unt of spiral stability over the lift-coefficient range 
and increasing the lift coefficient increased the spiral stability for 
both the sea— level and 30,000— foot-altitude conditions. 

The Air Force flying— qualities requirements do not call for spiral 
stability, but state that the spiral mode should not double amplitude in 
less than 1 seconds. It is therefore concluded that an airplane of this 
type will have satisfactory spiral stability for the conditions 
investigated. 


CONCLUDING EEMAEES 


The results of the investigation are summarized in the following 
paragraphs. Since small changes in some of the mass and aerodynamic 
characteristics may cause considerable differences in the lateral 
stability of an airplane, the results apply only to a 60° triangular- 
wing airplane in the condition for which the calculations were made. 


10 


NACA EM L9A04 


1. The lateral oscillatory stability was satisfactory over the lift- 
coefficient range for the airplane at sea level. 

2. Increasing altitude decreased the damping for all lift coef- 
ficients. The airplane was stable, at an altitude of 30,000 feet, but 
did not meet the U. S. Air Force specifications over the low— lift- 
coefficient range. 

3. The airplane was spirally stable for all conditions. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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TABU I 


E5ARACTKIBTTC8 Off A 60° IHXAKMAIWillO AIBEUK IBB) Hf IKKHStDCnHJ 

MS JWmttlrJAIIHAt^OSm^ BOmmABUS 

[b - 31.30 ft, B - 425. 00 b<l ft, V « 14114.0 It, y - 0°j i - 0.2l)0; 

Cv - 0, Cr - Ol Or. « -0.8 to 0.8; - 0.207 to 0. 

p r Ptfdl b 


Altitude 

Cl 

H 

C 

b 


*0 

% 

P 

Cy 

P(tail off) 

^{tau off) 

Ci 

Off) 

^tall off) 

c, 

^(tall off) 

^(tall off) 


m 

2.2J 

5.1 

-2,9 

13.83 

3.86 

9.48 

0.00238 

-0.11W 

-O.oss 

■-0.13 

0.005 

0.010 

-0.010 

Boo lerol 

.5 

n.6 

5.1 

6.5 

13.83 

3.86 

9.48 

0.00238 

-.200 

-.040 

-.10 

.030 

.018 

-.022 


1.00 

23.2 

?.i 

18.1 

13.83 

3.86 

9.48 

O.OOB38 

-.573 

-.100 

-.10 

.040 

-.172 

-.166 

■HAl 

31 

2.2 

5.1 

-e.9 

37.03 

3.86 

9.48 

0.000889 

-0.LU6 

-0.055 

-0.15 

0.005 

0.010 

-0.010 


.5 

11.6 

5.1 

6.5 

37.03 

3.86 

9.48 

.000889 

-.200 

-.040 

-.10 

.030 

.018 

-.022 

mm 

1.00 

23.2 

5.1 

18.1 

37.03 

3.86 

9.48 

.000889 

-.573 

-.100 



-.10 

.040 

-.172 

-.166 


Tfcd.1 ootiLrl tut Iona ana dotoraJjaed from the following a^patiooe: 


0; 

^-o) 

°i 

lr (Cnp-fr) 

° a P(Cnp-0) 

^^(Onp^O) 


ofr) t(-Vu „) 

C, r(*^ll crrr) + ^ "" t 0to “) (" Cn HaII off) 

✓ 

off) 

^(toll off) " S b (^V off) 


Cl “ Cv 

P b°ta 1 P (Cnp-0) 

Cs ■* O7 

* r tOtia ^(Cnp-O) 

C “ p totta ’ Cnp (Cap-0) 


-2 (f "t Blna ) 2 1 Cn p 

+ iS ljiBa )s 
+ s (b"t elna ) s 


- C„ 


total 


r (Cnp-0) 


- 8 H„ 
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TABIE H 


CTAEACTERISTICS OF A 60° TBIABGOIAE-WIRG AIRHAHE USED IN DSTESMTHING 
THE HSEIOD AND TIME TO DAMP TO OHE-HAUF AMFEETODE 


j^Tp "0| C Tr - 0; 7 - 0°] 


Altitude 

c L 


Ub 


(^ 2 

Kj 2 

/ 

KZ 2 

Exz 

p total 

\ 

. p total 

\ 

total 

C n 

^total 

c u 

r total 

n/pSV 

Sea level 

0.10 

•50 

1.00 

~e .9 

6.5 

18.1 

13.83 
— do- 
-do- 

0.0152 

— do — 
— do-- 

0.0916 

— do — 
— do — 

0.0154 

.0161 

.0225 

0.0914 

.0904 

.0842 

-0.00386 

.0086 

.0226 

- 0.1772 

-.1125 

-.1109 

-0.688 

-.715 

-.860 

0.0466 

.0396 

-.1437 

0 , 04 l 6 

.0515 

.0683 

-O.O59O 

-.0591 

-.2390 

0.821 

1,840 

2.590 

30,000 ft 

0.10 

.50 

1.00 

- 2.9 

6.5 

18.1 

37.03 

0.0152 

O.O916 

0.015li 

.0161 

.0225 


-0.00386 

.0086 

.0226 

- 0.1772 

-.1125 

-.1109 

-0.688 

-.715 

-.860 

0.0466 

.0396 

-.1437 

0.0416 

.0515 

.0683 

-O.O59O 

-.0591 

-.2390 

1.342 

3.001 

4.245 

— 0.0“ 
—do- 

--do — 
—do — 

— do— 




r t 


'■trovSi wa vove 





TABUS HI 

PERIOD AHD RECIPROCAL OP TIME TO DAMP TO ORE-HALF 
AMPLITUDE FOR A 60° TRIAIOILAR-WIIKr AIBPLAJIE 





Oscillatory mode 

Aperiodic (spiral) 
mode 

% 

(per deg)* 

(per deg) 

Period 

■ 91 

KpS 

1 

C l/2 

(l/cycles) 

1 

Tl/2 

(l/sec) 

O.OOO 83 

O.OOO 85 

2.97 

0.314 

0.935 

O.OO 79 O 

.000647 

.OOI 65 

4.66 

.361 

1.67 

.109 

.000912 

.OOO 953 

14.87 

> 0.66 

.389 

.389 

I .89 

8.06 


.OOO 83 

.OOO 85 

3.02 

. .141 

.426 

.00490 

,000647 

.OOI 65 

4.45 

.240 

1.06 

.0557 

.000912 

.000953 

A .76 

> 4.86 

.234 

.240 

1.11 

8.33 
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Figure 2.— System of axes and angular relationship in flight. Arrows 
indicate positive direction of angles, T] = 6 — 7 — e. 
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W/nc 

direction 


Wind into 

the page 


Wind 

direction 


Figure 3*— The stability-system of axes is defined, as an orthogonal 
system of axes having their origin, at the center of gravity and in 
vhich the Z-axis is in the plane of symmetry and perpendicular to the 
relative wind, the X-axis is in the plane of symmetry and perpendicular 
to the Z— axiSj and the Y-axis is perpendicular to the plane of symmetry. 
Aitovs indicate positive directions of moments and forces. 
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Figure 4.— Effect of lift coefficient on the neutral— lateral— oscillatory- 
stability boundary (R = 0) for tbe 60° triangular-wing airplane. Sea— 
level condition. 













Lift coefficient C< 
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